Non-local effects in the fermion Dynamical mean field framework. 
Application to the 2D Falicov-Kimball model. 
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We propose a new, controlled approximation scheme that explicitly includes the effects of non- 
local correlations on the D = oo solution. In contrast to usual D = oo, the selfenergy is selfconsis- 
tently coupled to two-particle correlation functions. The formalism is general, and is applied to the 
two-dimensional Falicov-Kimball model. Our approach possesses all the strengths of the large-D 
solution, and allows one to undertake a systematic study of the effects of inclusion of k-dependent 
effects on the D = oo picture. Results for the density of states p(u>), and the single particle spectral 
density for the 2D Falicov-Kimball model always yield positive definite p(w), and the spectral func- 
tion shows striking new features inaccessible in D = oo. Our results are in good agreement with the 
exact results known on the 2D Falikov-Kimball model. 
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PACS numbers: 71.27+a,71.28+d 

The effects of strong correlations on the properties of 
low-dimensional lattice fermion systems is still an open 
problem, inspite of many efforts spanning over thirty 
years. Recently, the development of a dynamical mean- 
field theory (DMFT), exact in D = oo has led to a major 
advancement in our understanding of the physics in the 
local limit DMFT is a nonperturbative scheme. 

It has provided a detailed picture of the Mott transi- 
tion in this limit, and has been fruitfully applied to a 
large class of models. Moreover, it has proved to be a 
rather good approximation to actual three-dimensional 
transition-metal and rare-earth compounds Inspite 
of its successes, DMFT has its shortcomings; the single 
particle selfenergy is k- independent, and so is not cou- 
pled to k-dependent collective excitations. This is clearly 
serious, e.g. near a magnetic phase transition where it 
cannot show any precursor effects. A consequence of the 
above is that the single particle spectral function can- 
not access such effects, and so the approach cannot be 
employed to describe e.g. angle-resolved photoemission 
experiments, or to study the changes in Fermi surface 
topologies driven by interactions in finite-!?. The above 
makes it imperative to develop such controlled extensions 
of DMFT as are able to rectify its unphysical aspects, 
while preserving its strengths. 

The DMFA maps the lattice system onto a selfconsis- 
tcntly embedded "impurity in a bath" problem that de- 
scribes the effects of the coupling to other lattice sites. As 
in a Weiss mean field theory, a self consistency condition 
is obtained by requiring the averaged Green's function of 
the bath Go to coincide with the local G at this impu- 
rity site. DMFA becomes exact in D — oo because one 
can show that the spatial correlations fall off at least as 
l/£>|i-j'l/2 with distance \i-j\ §. This "Weiss field" still 
has a nontrivial dynamics described by an effective local 
action Sdm f obtained by integrating out all sites exclud- 
ing the impurity site. Given Sdmf, one can compute 



the one particle Green's function on this site, Gi mp (u)), 
which is a functional of Go(w). The selfconsistency condi- 
tion that links Gi mp back to Go is solved iteratively. The 
main drawback of this method is that S is k- independent. 

In this letter, we present an extension of the DMFT 
that allows one to treat effects of non-local spatial cor- 
relations. To do so, we exploit the freedom involved in 
choosing the input bath propagator. We account for non- 
local correlations via the self consistent inclusion of two 
particle correlation function in the bath Green's func- 
tion (GF) Gq(uj). This is achieved by using the Spectral 
Density Approximation (SDA) jjj. Thus, our aim is to 
show that a proper combination of two methods -DMFT 
and SDA- is eminently suited to describe 1/D effects in 
lattice correlated fermionic systems. As far as applica- 
tion is concerned, we deal with the FKM as it was the 
first model to be solved exactly in D = oo. Attempts 
at a 1/D expansion for the FKM have also been made. 
Finally, there exist some exact results on the 2D FKM 
which are of great help to evaluate the quality of our ap- 
proximation, making it a first choice for testing methods 
in the context of extensions of the DMFT. 

We now present this new method. We begin with 
the formalism for the generic Hubbard model at zero- 
temperature. We then go to the FK to solve the full set 
of selfconsistent equations. We will focus on half-filling, 
the extension away from n = 1 is straightforward. The 
Hamiltonian is 



H 



<ij> : <7 i 



on a D— dimensional lattice. The Hubbard model is ob- 
tained when t-f = ti = t, and the FKM is given by 
t-f = t and ij = 0. In the case of half filling fi = U/2 
by particle-hole symmetry. These models have been ex- 
tensively studied using the DMFT. 

We get a bath propagator which includes 1/D effects 
using the SDA, first pioneered by Roth 0] to describe 
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magnetism in narrow bands. Its basic idea is to com- 
pute exactly the first few moments of the spectral den- 
sity to reconstruct the spectral function A(k,ui). With 
H as above, it is easy to compute the first four moments. 
They are given by the expectation value of the following 
?i-fold commutator 



=< [. 



-k,<r 



H],H],...H] > 



(1) 



Using usual scaling arguments Q, one sees that the 
higher n, the higher will be the powers of 1/D involved 
in the n-th moment. In the case of a generic Hub- 

(3) 

bard model, it turns out that the M^J contains all 



the o(l/ D) contributions 
o(w) + U 2n l( l ~ n i) B ^ with 
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It can be writen M^J 
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The SDA then permits one to write an explicit closed 
form expression for Aj(k, u>) = S[u> — ek + /i — So,-f(k, u/)], 
where £0 T C 4 ' w ) ^ s given ~~ for the generic Hubbard model 
-by 

E , T (k,^) = Un^U^l-n^uj+^-Uil-n^-B^]- 1 

This is the self energy in the SDA. It is interesting to 
notice that all the order 1/D effects enter only through 
-Bkcr, whereby the self energy acquires a non-trivial k- 
dependence. At this stage, this quantity is undetermined, 
and will have to be obtained self consistently. The SDA 
one particle Green's function can then be computed as 
G 0T (k,^) = [w-ek-Eo^k,^)]- 1 . 

We now consider the FKM, (tj = 0). In this case, the 
model has a local U(l) symmetry, [nn,H] = Vi. It 
thus follows from Elitzur's theorem that only the second 
correlator in B^ a is non zero. This leads to 



r H (I - ni )B kl = ~ Yl e- 4k -( R '- R ^[< n iinjl > -nf] 

<ij> 

= -iXf(k) (3) 



where XlC 4 ) is the order 1/D static structure factor of 
the down spin electrons. 

It is interesting to notice that the only order 1/D con- 
tribution to the single particle selfenergy is in fact 
The above equations are exact in the band as well as in 
the atomic limit, as can be easily checked. This has im- 
portant consequences, for it points to a way for further 
improvement of the formalism. It allows one to formulate 
a DMFT for the local part, more or less along same lines 
as in D = 00 0. However, in contrast to the regular 
DMFT where Eoo- = 0, Go contains information about 
non-trivial k-dependent features. Therefore, using this 



SDA Go a as the input one particle bath Green's function 
in a new DMFT scheme allows one to introduce in a sim- 
ple way the non-local spatial fluctuations at order 1/D 
in the FKM. 

We now come to the computation of the dynamical 
correlations. This results in a dynamically corrected self 
energy. To do this beyond the single-site level is a prob- 
lem that has not been attempted, and here we use the 
standard DMFT as an approximation to the 2D case. In 
the case of the half filled FKM, the method drawn from 
the lines sketched above can be developed the following 
way. 

As is known the equation-of-motion method 

(EOMM) gives the exact solution to the FKM in D = 00. 
We will use the EOMM to get the full (local) impurity 
Green's function which we denote by Gi mp ^(u>). 

The closed set of equations for the on-site Green's func- 
tion with a bath having Eo^k, u>) as a selfenergy is given 
in ref. @, replacing e k by e k + E , T (k,o;). 

The local Gj is computed to have the same form as that 
found in j2|, and the local self-energy is computed from 
Dyson's equation, E jmP:T (w) = G^ impA (uj) - G^ pA (u) 
where Go,imp,] {^) = [w + fi — A(w)] -1 , with 



A(u,)=£- 



'k 



M - £k - Eo,t(k,w) 



One then gets 



U 2 ni (l- ni ) 



+ M -C/(l-n i )-A(w) 



(4) 



(5) 



This local dynamical self energy is used to correct the 
local part of the bath self energy, to introduce dynamical 
effects in the k-dependent bath GF. This is done by 

E T (k, ui) = £ imP)T (u>) + £ , T (k,Lj) - ^ s o,t( k > w ) ( 6 ) 



and provides us with a dynamically corrected non-local 
Green's function G|(k, uj) = [u + — £k — E|(k, 
which is exact in both the atomic and the band limits. 
The density of states is then obtained by taking the imag- 
inary part of the k-summed G(k, ui). 

To complete the procedure, one should selfconsistently 
compute the order 1/D correlator which enters the bath 
Green's function as anounced earlier. In the case of the 
FKM, this means one has to compute the down spin 
static susceptibility. To begin with, we compute the 
vertex function exactly in D = 00. The full Xi(q) is 
then obtained from the Bethe-Salpeter equation in the 
framework of the EOMM @. The result is Xl(l) = 
ri|(l — n\/)/D{q), where 



^(q) = i-E 
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Gg + X ?(q) 
G?[x?(q)§7 + i] 



(7) 
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and x?(q) = (-1/N) Ek G t( k + q)G T (k), where the i/- 
sum is on the Matsubara indices (we take the zero tem- 
perature limit). In the above, the vertex corrections in 
the equation for the charge susceptibility of the down- 
spins electrons are approximated by their D — oo values. 
To compute the vertex function exactly to order l/D 
would be a formidable task, which has not been success- 
fully attempted even for the electron gas. 




(a) (b) 

FIG. 1. (a) The selfconsistency loop in usual DMFT. 
(b) Our method: k-dependent features are incorporated via 
an external input \l to the bath Green's function Go, which 
is computed self consistently as shown. 

To sum up, the equations (||)-(0) above form a com- 
plete set of selfconsistent equations for the FKM. These 
include explicit l/D effects through the k-dependence of 
£o,t(k, oj). To solve these equations, we start with an ar- 
bitrary XJ.(*l) as an input in our initial Got- This is used 
to compute Eo,f(k, to), and the full Gi mv ,\ (u>) from the 
large- D solution along the lines of 0. At this step, we 
have a DMF approximation to the FKM that is equiva- 
lent to a selfconsistent loop, with an external input quan- 
tity, the structure factor. The second step of our method 
is to compute Xl(q)- F° r this purpose, we use the DMF 
approximation to the structure factor equation (Q), plug- 
ging the obtained Gj (k, to) along with the full local ver- 
tex function in the Bethe-Salpeter equation to compute 
Xx(q)- This is fed back into the equation for Eo,j and 
the whole process is iterated to convergence. The differ- 
ence between the regular self consistent DMFT and our 
method is illustrated in fig.l. 

We implemented numerically this new method to study 
the Falicov-Kimball model on a 2D square lattice at 
T = 0. The up spin electron dispersion relation is 
thus ek = — 2t(cosk x + cosfc y ), and the corresponding 
density of states has a van Hove singularity at w = 0. 
The k sums are performed on a L x L mesh, while 
the w-integrals are computed using standard integra- 
tion procedures, using a small broadening 77 to repre- 
sent 8 functions. Convergence is obtained when \ 2 — 
J2uj \Gnew{u) — G id{u)\ 2 < £ where e is a small thresh- 
old. At each step of the iterative process, a check on 
the accuracy of the computation is provided by the Lut- 
tinger sum-rule ^ k J_ duM CT (k, lu) — n a , which must 



be satisfied to high accuracy. 

We now describe the results of our computation. 
Working with a lattice size L = 64, and a threshold 
e = 10 -9 , we observe a very fast convergence of the dou- 
ble selfconsistent procedure, that needs less than 10 steps 
to reach the convergence threshold for Gf (u>). The Lut- 
tinger sum rule is satisfied with an error of less than 10~ 3 . 




FIG. 2. DOS for the 2D Falicov-Kimball model for sev- 
eral values of U/t at half-filling, with our method. It is always 
positive definite. The inset shows the behavior of the U/t — 2 
gap at lu = when changing the artificial broadening 77 of the 
S functions, rj/t = 5.1CF 2 thick line, Tj/t = 5.10 -3 thin line 
and r//t = 2.10 -4 dashed line. 

Fig. 2 shows the density of states (DOS) p(u>) obtained 
by our method for several values of U /t. Firstly, notice 
that this function is always positive definite, in contrast 
to results of rcf. flqj. This is a positive feature of the 
present method. We should stress that we are interested 
in the effect of coupling of self-energy to l/D fluctua- 
tions, which, for the FKM, are completely static. The 
effect of the interaction is very clear. Starting from the 
2D free fermions value at U/t — 0, it opens a gap at 
lu = 0. The critical (U/t) c for opening up a gap is very 
small. As a matter of fact, working with a L x L lattice we 
cannot resolve energy details smaller than 2-kvf/U ~ 0.2 
here. However, playing with the parameters of the prob- 
lem, we see a real gap with U /t as low as 0.1. This result 
agrees with the exact result on the 2D FKM at weak cou- 
pling ly}, where the gap and the related checkerboard 
order arises from the nested FS at n = 1. The inset of 
fig. 2 shows the behavior of the DOS at low frequency 
for U/t = 2, obtained by changing the broadening of the 
delta functions to get a better precision. The features 
observed at higher frequencies are related to the unper- 
turbed 2D band structure and are absent in infinite D. 

The computed Xi(l)j which contains information 
about long-range order (LRO), shows interesting behav- 
ior; along the zone-diagonal, at the point (k = [tt, 7t]), it is 
a constant, scaling approximately with lattice size. This 
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is the expected behavior in the broken-symmetry phase 
]l2fl , and implies that the n = 1 ground state is charac- 
terized by checkerboard LRO of the down-spin electrons 
||. In ref. [O, the behavior of Xi(q) can be studied only 
in ID, because of finite size restrictions. In contrast, our 
approach yields unambiguous conclusions in the thermo- 
dynamic limit in 2D. Again, this is in clear agreement 
with the exact result [jll| . 




FIG. 3. The spectral function ^4 T ( k ,^) for tne 2D Fau - 
cov-Kimball model at half-filling, for U/t = 4, and k along 
the zone diagonal, k = (mr/8,mr/8), n going from to 8. 

The greatest advantage of our method, however, is 
that it allows a detailed study of k-dependent (at the 
order of 1/D) effects in the one particle spectral func- 
tion A(k, oj) = — ImG(k, oj)/it. In fig. 3, we show A(k, oj) 
for k along the zone-diagonal from k = to k = (it, tt) 
for a representative value of U = 4t. We notice that 
A(k, oj) satisfies all the symmetry properties consistent 
with those expected from pure nearest neighbor hopping 
on bipartite lattices, and with particle-hole symmetry: 
v4(ko,w) = A(ko,— oj) for ko = [tt/2,tt/2], while along 
the zone diagonal A(ko +5, id) = A(ko — 8, —oj), V<5 along 
the zone diagonal. We show only the A(k, oj) for oj < 0; 
the oj > part can be obtained from the symmetry prop- 
erty mentioned above. We have checked that our calcula- 
tions are fully consistent with the above property. More- 
over, we see that since Q = [k, tt] is a nesting vector in 
2D, the above property of A(k, oj) with 5 = [tt/2,tt/2] 
leads to a natural explanation of the "shadow-band" fea- 
tures observable in ARPES. It would be very interesting 
to see whether they survive away from n = 1 . 

It is instructive to compare our results to those ob- 
tained by Velicky et. al [ flO| in their pioneering paper on 
CPA. CPA is equivalent to the exact D — oo result for the 
FKM, and so a comparison with |l(J allows us to study 
the effects of 1/D effects on A(k, oj). A comparison of our 
data with results for identical parameters from p0[ shows 



that the coupling to 1/D fluctuations induces new fea- 
tures in A(/k, oj) compared to those observed in DMFT. 
In D = oo, A(k, oj) has a two-peak structure with the 
dispersion controlled solely by the free bandstructure. In 
our approach, this is modified because of the extra k- 
dependence coming through the -Bki- This shows that 
a formalism capable of explicitly treating intersite cor- 
relations permits one to access k-dependent features in 
A(k, oj), so we suggest that this method can be fruitfully 
applied to compute ARPES lineshapes in correlated sys- 
tems. A more detailed application to spectroscopy would 
require the use of the actual bandstructure DOS, and is 
left for future work. 

In conclusion, we have developed a simple, physically 
appealing way to study the effects of non-local spatial 
fluctuations on the single particle spectral properties. We 
have applied the formalism to compute the single parti- 
cle spectral function and the DOS for the 2D FKM and 
have obtained results in good accordance with what is 
known from the exact solution. Extensions of the work 
to look at the metallic phase in the FKM off n = 1, as 
well as for the Hubbard model, are being studied and will 
be reported separately. 
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